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The perfect-nesting instability towards antiferromag- 
netism of the Hubbard model is suppressed by next-nearest 
neighbor hopping t' . The asymptotic behavior of the crit- 
ical coupling Uc{t') at small t' is calculated in dimensions 
d = 2, 3, cxD using Hartree theory; this yields the exact result 
at least in d > 2. The order of the transition is also deter- 
mined. A region of stability of a metallic antiferromagnetic 
phase in d = 3 is identified. 



I. INTRODUCTION 

Microscopic investigations of strongly correlated elec- 
tron systems are usually based on lattice models, the sim- 
plest of which is the one-band Hubbard model. In this 
model two particularly important features of real materi- 
als, the structure and the dimension of the underlying lat- 
tice, enter only via the kinetic energy. In the past, most 
model calculations started from a given lattice (mostly 
hypercubic) and, for simplicity, assumed the hopping to 
be restricted to nearest neighbors (NN). In this case lat- 
tices are either bipartite or not. In dimensions d > 1 
and at half filling bipartite lattices generically lead to an 
antiferromagnetic instability in the ground state for any 
U > due to Fermi surface nesting. 

It has recently become clear, however, that pure NN- 
hopping with amplitude t cannot account for some impor- 
tant features of correlated systems and that next-nearest 
neighbor (NNN) hopping with amplitude t' has to be in- 
cluded for both qualitative and quantitative reasons. In- 
deed, realistic band structures, photoemission data and 
neutron-scattering measurements of high-Tc and related 
materials usually cannot be fit by an epjergy dispersion 
that originates solely from NN-hoppingJcl; estimates for 
\t'/t\ range from 0.15 to 0.5. A particularly clear-cut 
example is the antiferromagnetic insulator Sr2Cu02Cl2 
where the inclusion of a t'-term is essential to understand 
the ARPES datacl. The hole motion in|-|the f-J-model is 
also greatly affected by NNN-hoppingJ. Furthermore, 
it became clear only most recently that the stability of 
metallic ferromagnetism in the one-band Hubbard model 
strongly dependsr-an the presence of a i'-term in all di- 
mensions d > IB^Ilil. 

On bipartite lattices in d > 2 the t'-term has pro- 
found consequences: the frustration introduced by t' sup- 
presses the perfect-nesting instability towards antiferro- 
magnetism at small U . Up to now_this effect was mainly 
studied in d = 2. Lin and Hirschlij calculated magnetic 
properties within Hartree theory and by exact diagonal- 



ization of small systems. Tremblay et al.0 studied both 
magnetic and pairing correlations by quantum Monte- 
Carlo simulations and a tmo-particle self-consistent ap- 
proach. Duffy and MoreoEj investigated the renormaL 
ization of t' by the interaction U, and presented resultscil 
suggesting the existence of an itinerant antiferromagnetic 
ground state at half-filling. Metallic antiferromagjietic 
phases areJuiown to exist, for example, in V2-x03liJ and 
NiS2-xSeJia. 

In this paper we wish to obtain analytic insight into 
the suppression of the antiferromagnetic nesting insta- 
bility at small U and half filling due to the frustration 
introduced by t'. To this end we solve the t-t' Hubbard 
model 

^ = - E * (^Is- + ^■'^•) + E + '^■c-) 

<ij> (ij) 

within Hartree theory at T = and thereby calculate 
Uc{t'), the critical value of U for the onset of antiferro- 
magnetism. FaCjSmall t' the results for Uc become asymp- 
totically exactO, at least in dimensions d > 2. 

We will calculate Uc{t') explicitly in d — 2,3,oo, both 
numerically and analytically and will also determine the 
order of the transition into the antiferromagnetic state. 

II. HARTREE EQUATIONS 

The Hubbard interaction is purely local and there- 
fore does not include any direct exchange. Conse- 
quently, there is no Fock term and Hartree-Fock the- 
ory reduces to Hartree theory. Within that approach 
spontaneous symmetry breaking with an antiferromag- 
netic order parameter (Neel order) is introduced by the 
ansatz (riio-) = ^(1 + 1)'*' A) with a — ±1, where 
A is the staggered magnetization. We can then de- 
compose the electron number operator into its expecta- 
tion value and a fluctuating part : : according to 
riio- —: Uia- : -f (n^o-). Neglecting terms quadratic in the 
fluctuations, the Hartree Hamiltonian takes the form 
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where the summation is restricted to half of the Brillouin 
zone. Here we have introduced the number of sites N and 
the nesting vector Q = (tt, . . . , tt). The dispersion of the 
free electrons on the lattice is given by 



Ek = — 2t cos ki + At' cos ki cos kj 



(3) 



For t' ^ the perfect nesting property e(k) = — e(k + Q) 
of the free dispersion no longer holds. In the following we 
shall use the notation k' = k + Q. After diagonalization 
of the above quadratic form the mean-field Hamiltonian 
can be written as 



UNA^ 



(4) 



with the single particle energies 



ek=2 



[/2A2 

(ek-Ck'^ 



(5) 



Since ek does not fulfill the perfect nesting condition, one 
has to adjust the chemical potential as a function of (3, A, 
U and t' to ensure the correct filling n — I. In addition, 
the value of the order parameter A is determined self- 
consistently by the condition that the free energy 



F[P, A, U] = -p-^ In Tr exp(-/37^) + ^iV, 



(6) 



be minimal with respect to A. This leads to two self- 
consistency equations: 
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(27r)" 
d'^k 



f(/3, ck - m) 



(7) 

(8) 



where f(/3, ek — /i) is the Fermi function. The trivial solu- 
tion A = of the minimum condition has been excluded 
in (||). In the following, (^) and will be analyzed for 
dimensions d = 2, 3 and cxd at T = 0. 



III. RESULTS AND DISCUSSION 
A. d = 2 

We set t = 1. The presence of i' 7^ prevents an 
analytic evaluation of the integrals in (|^) and (||). We 
therefore computed them numerically by k-summation 
on reciprocal lattices of sizes up to 8000 x 8000. Intro- 
ducing the notation 



7 := 



UA 

2 



(9) 



for the energy gap of the Hartree dispersion , we notice 
that both the integrals in (0) and (||) only depend on 7. 

We find that for any t' ^ 0, the integral in (^) has a 
finite limit for 7 — > 0. This implies that below a critical 
value t/c of the interaction a stable solution with A > 
no longer exists—For d = 2 this fact was already noted by 
Lin and Hirscht3. The dependence of Uc on t' is shown 
in fig. |. 
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FIG. 1. Critical interaction Uc vs. next-nearest neighbor 
hopping matrixelement t' in d = 2. Dashed line: asymptotic 
behaviour, (pj|). 

For t' 0, Uc{t') decreases very slowly. An analytic 
description of the behavior of Uc{t') may be obtained 
by noting that the divergence of the integral in (^) for 
t' = iSrCut off just at the value of 7 where the indirect 
bandgaplia 




FIG. 2. Hartree dispersion ek in d = 2 for t' — 0.4 and 
7 — 1.0 with an indirect bandgap between (7r/2,7r/2) and 
(0,^). 
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in the Hartree dispersion (||) closes, namely at 

7* = 2t'. (10) 

The asymptotic behaviour of Uc{t') for small t' may thus 
be found by setting t' = in (^) and using 7* as a cut-off. 
This leads to 



[Ucit')]-^ ^ ci / deK 
Jo 



where K is the complete elli ptic integral and ci = 




l/(27r^). The approximation (|11 
t' ~ 1, as-can be seen from fig 
reduces tolZl 



is very good up to 
For sniaU t', (IF 



mt')]- 



Cl 



(Int') +C2 Ini' 



(12) 



where C2 ~ —0.15. The quadratic dependence of Uc on 
\nt' is characteristic for d = 2, where the logarithmic 
divergence of the density of states (DOS) for i' = at 
the Fermi level introduces an additional factor Int' . 

The behaviour of the order parameter A at U Uc 
can be characterized as follows: For t' below a threshold 
value ip ~ 0.4 the order parameter A jumps to discon- 
tinuously when the interaction is reduced below Uc{t'), 
indicating a first order phase transition. For U > Uc, 
the order parameter takes exactly the same value as for 
t' = 0. For t' > and U ^ Uc the order parameter goes 
to zero continuously, but extremely fast. 



B. d = 3 

We evaluate the Hartree equations numerically at zero 
temperature, using lattice sizes up to 400 x 400 x 400. 
As before we notice that the divergence of the integral in 
(^ for 7 ^ is cut off for any finite t'. The resulting 
critical interaction Uc{t') is shown in fig. O. 




FIG. 3. Critical interaction Uc vs. t' in d — 3. Dashed 
line: asymptotic behavior, ([l^). 



The indirect bandgap again closes at 7* = 2t' , but in 
this case the order parameter goes to zero continuously 
as U Uc, indicating a second-order phase transition. 
The logarithmic divergence of the integral in (||) is, cut 
off at 7*, which leads to the asymptotic cxpressionO 



[Uc{t')]- 



'C3 Int' 



C4, 



(13) 



where C4 « 0.25 is a fit parameter and C3 = N{0) « 0.14, 
with A^(0) as the DOS at the Fermi level in d = 3 for 
t' = 0. Eq. (|l|) represents a good fit for t' < 0.4. 



C. d = 00 

To analyze the model (|l|) on a hypercubic lattice in the 
limit of infinite dimensions, one has to ensure thalj— thsj 
hopping remains finite; this is achieved by the scalinglj'E^I 



t* 



2d 



t' 



t' 



y/2d{d~ 1) 



(14) 



for the NN- and NNN-terms, respectively. In the ijollow- 
ing we set t* = 1. As shown by Miiller-HartmannE] the 
dispersion of the free electrons then reduces to 



t' 



(15) 



Here e is the kinetic energy|-without the t'*-term, which 
has a Gaussian distributionEj 



Nie) 



1 



(16) 



As a result, we can write the Hartree dispersion in the 
form 



-V2t'* (e2 - 1) 



2eA 1 



C/2A2 



(17) 



as before, it only depends on t'* and 7 = UA/2. We 
note that due to the exponential band-tails of the hyper- 
cubic DOS there is no true gap between the two Hartree 
bands for any value of 7. The Hartree equations can 
be written as one-dimensional integrals and are readily 
solved numerically. Again we find a critical interaction 
Uc for the stability of long-range antiferromagnetic or- 
der at T = 0. A plot of Ucit'*) is shown in fig. |. Due 
to the absence of a band gap, the asymptotic behavior 
of Uc for small t'* is qualitatively different from that in 
d = 2, 3. In the present case, the cut-off for the integral 
in (^) is provided by a small filling in the upper one of 
the two Hartree bands (|l^) which vanishes exponentially 
for t' ^ i.e. 7* w exp(— ^^,lyi ). This leads to an 
asymptotic behaviouJlZI of the form 



Uc = c{t'*f 



(18) 
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where c « \J2t^. Hence Uc decreases much faster than 
in d = 2, 3. As in d = 3 the order parameter A goes 
to zero continuously for U ^ Uc- A second-order phase 
transition is therefore seen to be generic for d > 2. 
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FIG. 4. Critical interaction versus t'* in d = oo. Dashed 
line: asymptotic behavior, 

We note that in d = 3 , the continuously vanishing 
order parameter implies a (small) range of U values where 
the band gap is zero but A is still finite. This implies the 
existence of a metallic antiferromagnetic phase for i' 7^ 
as shown in fig. 5. 
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Note added in proof: 

In d = 2, if the i'-hopping occurs only along one 
of the diagonals of the squares, the dispersion reads 
Ek = — 2t (cos kj; + cos ky) + 2t' cos{kx + ky). The corre- 
sponding t— i'-Hubbard model may be viewed as a simple 
dimer model of the two-dimensional organic supercon- 
ductor K - (BEDT - TTF)2X which was investigated by 
Kino and Fukuyama [J. Phys. Soc. Jpn. 65 (1996) 2158] 
within Ha rtree- Fock theory. In contrast to the results 
of section HI A , one then obtains a region of stability 
of an AFM phase which is qualitatively similar to the 
result for the regular t — i'-model in d = 3 shown in 
fig. 5. The discontinuous transition from the paramag- 
netic metal into the antiferromagnetic insulator obtained 
by Kino and Fukuyama is found to occur for \t'\ > 0.6 
when this AFM phase is unstable. 

We thank R. McKenzie for drawing our attention to this 
variant of the t — i'-Hubbard model, and to H. Kino and 
H. Fukuyama for correspondence. 
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FIG. 5. Region of stability of a metallic antiferromagnetic 
phase (AFM) in d = 3 in the U-t' phase diagram. Full line: 
Uc{t'), dashed line: value of U where the band gap closes. 
Outside the AFM region the system is either an antiferro- 
magnetic insulator (AFM) or a paramagnetic metal (PM). 

In d = cxD, on the other hand, there is never a true 
bandgap anyway, as described above. Therefore, for 
t' ^ and U > Uc we always find metallic antiferro- 
magnetism. The existence n£.such a phase was already 
discussed by Georges et al.EJ. Quite generally, in di- 
mensions 2 < d < 00 a next-neighbour hopping term 
t' is found to be very effective in suppressing the perfect 
nesting instability towards an antiferromagnetic ground 
state. 
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